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Abstract 



A brief survey of the theoretical, numerical and experimental studies of the random field Ising model 
during last three decades is given. Nature of the phase transition in the three-dimensional RFIM with 
Gaussian random fields is discussed. Using simple scaling arguments it is shown that if the strength of 
the random fields is not too small (bigger than a certain threshold value) the finite temperature phase 
transition in this system is equivalent to the low-temperature order- disorder transition which takes place 
at variations of the strength of the random fields. Detailed study of the zero-temperature phase transition 
in terms of simple probabilistic arguments and modified mean-field approach (which take into account 
nearest-neighbors spin-spin correlations) is given. It is shown that if all thermally activated processes are 
suppressed the ferromagnetic order parameter m(h) as the function of the strength h of the random fields 
becomes history dependent. In particular, the behavior of the magnetization curves m(h) for increasing 
and for decreasing h reveals the hysteresis loop. 
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1 Introduction 

The thermodynamic properties of the random field Ising model (RFIM) remains the controversial issue during 
more that 30 years now (for reviews see e.g. [1,2]). Traditionally there are exist three somewhat independent 
mainstreams of reseach: theoretical, numerical and experimental. Inside each of them, as well as between 
any of them, one finds broad divergences in the obtained results and in their interpretations. The main issue 
of this controversy is the structure of the phase diagram of the system: what kind of phase transitions and 
which thermodynamic phases are present there? 

The random field model has been proposed originally by Larkin [3] to model the defect pining of vortices 
in sumerconductors. The most simple version of the this model for systems with discrete Ising symmetry, the 
random field Ising model (RFIM), can be described in terms of the Hamiltonian 



where the Ising spins {<Ji = ±1} are placed in the vertices of a D-dimensional lattice with the ferromagnetic 
interactions between the nearest neighbors, and quenched random fields {hi} are described either by the 
symmetric Gaussian distribution with {hihj} = 5ijh 2 , or by the binary distribution in which the fields are 
taking two values ±/i with independent probabilities at each site. 

1.1 Theory 

The first round of controversies in the studies of RFIM was about its lower critical dimension. According 
to simple energy ballace physical arguments by Imry and Ma [4] it should be expected that the dimensions 
D c above which the ferromagnetic ground state is stable at low temperatures (it is called the lower critical 
dimension) must be D c = 2 (unlike the pure Ising systems where D c — 1). Later, the existence of the 
ferromagnetic long range order in the 3D RFIM has been confirmed by a rigorous proof [5]. On the other hand, 
a perturbative renormalization group study of the phase transition demonstrates the so called phenomenon 
of a dimensional reduction, such that the critical exponents of the random field system in the dimension D 
appear to be the same as those of the ferromagnetic system without random fields in the dimension d=D-2 
[6], which implies that the lower critical dimension of the random field Ising model must be equal to 3, 
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in contradiction with the rigorous results. Actually, the procedure of summation of the leading large scale 
divergences assumes that the Hamiltonian has only one minimum. However, one can easily see that as soon 
as the temperature is close enough to the putative critical point (as well as in the whole low temperature 
region), there are local values of the magnetic fields for which the free energy has more than one minimum 
[7,8], and therefore the dimensional reduction is not grounded. Thus, the above arguments have settled the 
controversy about the lower critical dimension of the random field Ising model in favor of the value D c = 2. 

Phase diagram in the (h,T) plane for the RF Ising system has been first derived in terms of simple 
mean-field theory for the model with infinite-range ferromagnetic interactions (which corresponds to infinite 
space dimensionality) with the Gaussian distribution of random fields [9]. In this system the ferromagnetic 
phase is separated from the paramagnetic phase by a line T(h) of the second order phase transition which 
monotonically decreases with h from the point T c (0) of the pure system, down to the ending point T = at 
h = h c . On the other hand, for the bimodal distribution [10] one finds that the line T c {h) contains a tricritical 
point (T*,/i»), such that at h > /i* the phase transition becomes the first order. Later on a simple scaling 
description of the RFIM phase transition in finite dimensions has been developed assuming that it is of the 
second order, and that it is controlled by the zero-temperature fixed point [11]. 

As another extreme case, instead of the Ising spins, one can study the m-component vector spin system 
(the Ising model corresponds to m = 2) with Gaussian random fields. It can be shown [12] that in terms 
of the replica field theory in the limit m — > oo the statistics of interfaces in this system is described by the 
solution with broken replica symmetry (similar to that in spin glasses [13]) which is consistent with the low 
critical dimension di = 2. Later on in the same limit m — > oo it has been demonstrated [14] the instability 
of the replica symmetric state and the onset of the replica symmetry breaking scenario. In terms of this 
approach it has been argued that in the phase diagram of the three-dimensional RFIM the the paramagnetic 
and ferromagnetic phases are separated by the glassy phase in which the replica symmetry is broken [15]. 

An independent study of the <j) A replica field theory (which in the critical region is believed to describe 
RFIM) in terms of the Legendre transforms technique and the virial-likc relations also shows the presence of 
the intermediate RSB glassy phase below dimensions D = 6 [16]. Further renormalization-group (RG) study 
of this field theory in dimensions D < 6 revealed no stable fixed points signalling the nonset of the replica 
symmetry breaking [17]. 

Another line of theoretical studies concerns the non-perturbative thermodynamic states which are ap- 
parently missing in the usual RG calculations and which may appear to be quite relevant for the nature of 
the phase transition in RFIM. The importance of the non-perturbative phenomena has been noted already 
long time ago [7, 18, 19]. In physical terms, the non-perturbative configurations are rare spatial regions with 
"flipped" (opposite to the background) magnetizations, and in terms of the replica field theory they are de- 
scribed by localized in space finite energy instanton configurations in which the replica symmetry is broken. 
In the recent investigations a systematic approach for the calculation of the non-perturbative contributions 
has been developed [21] and it has been shown that from the point of view of their relevance the dimension 
D = 3 turns out to be marginal: formally these degrees of freedom produce finite (non-analytic in the strength 
of the field h) contribution only in dimensions D < 3 [22] . 

Finally, the construction of an alternative approach, the functional renormalization group, which is sup- 
posed to take into account all non-perturbative degrees of freedom for the whole class of the random-field 
O(N) spin model is under way these very days [23] 

1.2 Numerics 

Numerical investigations of the phase transition in RFIM reveal broad controversy in the obtained results 
and in particular in their interpretation. The main point of the disagreements is about the nature if the phase 
transition: is it continuous, or is it of first order with finite jump of the order parameter. 

Originally, the analysis of the high-temperature series expansion for the RFIM with the Gaussian distri- 
bution of random fields [24] indicated on the existence of a fluctuation-driven first order phase transition at 
sufficiently strong disorder below four dimensions (which suggests the existence of a tricritical point in the 
three-dimensional system). However more recent study indicated on the continuous transition, at leas for 
weak disorder [25]. 

First extensive Monte Carlo (MC) numerical simulations for the 3D RFIM have been carried out by Young 
and Naucnbcrg [26], who reported that the transition may be first order. On the other hand, the numerical 
study by Ogielski and Huse [27], and by Ogielski [28] shows that the transition may be second order. More 
recent Monte Carlo simulations [29] can also been interpreted as describing a continuous transition but with 
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a finite jump in the magnetization. However, later extensive MC calculations [30] for the energy and the 
magnetization distributions suggest a first order transition. 

The zero-temperature studies play an important role in understanding the nature of the phase transi- 
tion in the RFIM. According to the zero temperature fixed point hypothesis, the transitions which take 
place at T = and at T ^ are in the same universality class. First zero-temperature Migdal-Kadanoff 
renormalization-group calculations indicated a continuous transition which is characterized by very small 
value of the magnetization critical exponent (/3 ~ 0.02) [31]. Similar results were obtained in the finite tem- 
perature renormalization group study [32] . On the other hand, some of more recent zero-temperature finite 
size scaling calculations [33] have demonstrated the existence of a first order phase transition, although most 
of the them [34, 35] are in favor of the continuous transition. Moreover, in one of the most recent calculations 
it has been demonstrated that for a given realization of disorder the ground states (T = 0) and thermal states 
(T =/= 0) near the critical line are strongly correlated, which may be interpreted as the concrete manifestation 
of the zero temperature fixed point scenario [36]. The most precise results obtained by A.A.Middleton and 
D.S.Fisher [35] are clearly consistent with the existence of a single phase transition (no intermediate spin-glass 
phase has been observed) and it is characterized by very small value of the magnetization critical exponent 
= 0.017. 

Others numerical studies of the phase diagram of the RFIM in the plane (T, h) [37, 38] , also find no glassy 
phase (with or without replica symmetry breaking) which would separate paramagnetic and ferromagnetic 
phases. On the other hand, for the bimodal random field distributions a discontinuity in the magnetization 
at the phase transition line is observed [38] . 

Finally, the numerical simulations of the diluted antifcrromagnct in a field (DAFF) designed to imitate 
the random field magnetic systems which are studied experimentally (see below), shows that in equilibrium, 
the transition to the LRO state may be first order [39,40], and moreover it is claimed that at sufficiently 
strong dilution there is first order phase transition to a spin glass state [41] 

To conclude this brief survey of numerical studies, it should be noted that all numerical simulations in 
RFIM are impeded by the dramatic slowing down when approaching a putative phase transition, which very 
often makes very difficult to interpret the obtained results. Besides, the delicate point in all these studies is 
that the order parameter critical exponent (3 turns out to be very small, so it is very difficult to distinguish 
the second order transition (for which the order parameter continuously goes to zero at T c ) from the first 
order (when the order parameter has a finite jump at T c ). On top of that, recent numerical study of the 
critical region of the 3D RFIM indicates that the situation could be even more sophisticated due to absence 
of self-averaging of the correlation length, of the specific heat, and maybe of some others thermodynamical 
quantities [42]. 

1.3 Experiments 

From the point of view of experimental investigations an important theoretical breakthrough occurred when 
it has been realized that RFIM can be generated in diluted Ising antiferromagnets by application of a uniform 
external magnetic field [43]. This has opened the way for experimental studies of RFIM, and in particular 
for the investigation of the phase transition. 

In the very first experimental studies of the RFIM performed in the diluted antiferromagnets Co x Zni- x F 2 
it has been claimed that even the smallest magnetic fields destroy the long-range order at all temperatures 
in three dimensions [44]. 

First convincing experimental demonstrations of the existence of the phase transition in the 3D random 
field antiferromagnet Fe x Zn\- x Fi and its critical properties (assuming that it is of second order) has been 
reported in [45]. The extensive study of the critical properties of the system Feo.eZno.4F2 for temperatures 
T > T c (h) has been done in [46]. 

In fact, the behavior observed near the critical point depends on the field and temperature procedures 
used in the measurements. In the field cooled (FC) samples the transition appears to be rounded, and at low 
temperatures the system freezes in the metastable domain state. However, if the sample is cooled in the zero 
field (ZFC), the transition appears to be much sharper, and rounded only by slow dynamic effects [2]. 

On the other hand, in another experimental study of the 3D RFIM (for weakly anisotropic Mno.5Zno.5F2) 
[47] it has been observed that when approaching the phase transition the correlation length reaches a finite 
size, at which point it was assumed that further approach towards a putative second order transition is 
interrumpted by occurrence of a first order transition. 

The most slippery issue of the experimental studies of the phase transition in the RFIM is the measurement 
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of the static critical behavior of the staggered magnetization. Similarly to the numerical studies, the fact 
that critical exponent (3 is experimentally very small (e.g. according to [48] (3 ~ 1/8 or less, according to [49] 
/3 ~ 0.16) makes it very difficult to work out any definite conclusion about the nature of the observed phase 
transition. 

In the studies of the diluted antiferromagnet Mno.75Zno.25F2 the ficld-coolcd data were interpreted as the 
evidence for an equilibrium second-order transition sXTpc, while the zero-field cooled data were interpreted 
as non-equilibrium trompe I'oeil transition in which the long-range order diminishes continuously at Tzfc > 
Tpci while the correlation length only reaches a finite maximum value at a temperature between Tzfc and 
Tpc [50]. Similar behavior was observed in Feo.5Zno.5F2 [51]. On the other hand, all these results can be 
explained as the mean-field first-order transition broadened by the cluster- flipping mechanism [52]. 

In the study of the random field antiferromagnet Feo.5Zno.5F2 a notable broadening of the phase transition 
region has been observed, so that the data may be described by the Gaussian distribution of effective transition 
temperatures with the width which scales with the strength of the applied field as h 2 [53]. Besides, the 
apparent critical behavior observed in these measuremets represented a continuous evolution from metastable 
behavior towards equilibrium behavior. Moreover, according to the point of view of Birgencau [54], so far no 
true equilibrium phase transition has been observed in RFIM systems. Nevertheless, in the recent study of 
high-magnetic concentration Ising antiferromagnet F e . 93 Zn n . 07 F2 (which does not exhibit the severe critical 
scattering hysteresis [55]) it is claimed that in terms of general scaling assumptions the critical scattering 
analysis allows to reveal equilibrium critical behavior and to obtain the corresponding critical exponents [56] 
(in particular the specific heat critical behavior is well fitted by the logarithmic singularity which corresponds 
to a = [57]). 



In conclusion of this brief historic review one can note that the situation with the phase transition in the 
three-dimensional RFIM remains far from being clear. Theorists are mostly interested in their own problems 
(like 1/m expansion, or RSB in dimensions close to 6, or two-loop approximation for O(N) sigma model 
with the number of spin components close to iV c ~ 2.835), which are rather far from the original system, 
while people doing numerics and in particular experiments are facing realities full of numerous secondary and 
accessory effects which very often overshadow the main physical phenomena. 

It is evident that one can not hope to find the exact solution for this problem, and therefore some kind 
of simplifications and approximations are unavoidable. In this paper, which is mostly devoted to the low- 
temperature part of the phase diagram, it will be supposed that: (a) the finite temperature thermodynamical 
properties of the three-dimensional RFIM can be described in terms of the continuous Ginzburg-Landau </> 4 
field theory; (b) we can neglect the presence of randomness in the spin-spin interactions; and (c) in the 
study of the zero-temperature phase transition (at variation of the strength of the random fields) all thermally 
activated processes are suppressed. 

The above reservations are rather essential. 

(a) Ginzburg-Landau field theory. We all believe in the universality, and therefore if we admit 
that the phase transition in a given system is of the second order, then it is generally accepted that its 
critical properties can be described in terms of the corresponding Ginzburg-Landau Hamiltonian (with proper 
symmetry properties). However, if the system is considered at a finite distance from the putative critical 
point (which is the case in the present paper) where the correlation length is finite, then the relevance of 
the continuous limit Ginzburg-Landau Hamiltonian for the original lattice Ising system becomes much less 
evident. 

(b) Randomness in the spin-spin interactions. Apart from the numerical investigations of DAFF 
[39,41], the effects of randomness in the spin-spin interactions is practically never taken into account in 
theoretical studies (there are enough headache due to random fields...). However, it is well known that the 
presence of this type of randomness (without random fields) may change the critical properties of a magnetic 
system, and in particular it definitely does in the three-dimensional Ising model. Moreover, in some statistical 
systems (e.g. in the Potts model) the presence of such disorder may turn the first-order transition (of the pure 
system) into the second order one [58]. Therefore, it we admit for the moment that the phase transition in the 
"pure" (without randomness in spin-spin interactions) 3D RFIM is of the first order, it does not guarantee 
that in the corresponding experimental realizations of this system (where the randomness in the spin-spin 
interactions is inevitable) the phase transition would not turn into the second order one. 

(c) Thermal activations. In the studies of the zero-temperature properties of the system as the function 
of the effective strength of the random fields h, one can consider two types of procedures. In the first one it 
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is assumed that at any given h the system is at the ground state (which, at a given realization of the random 
fields, is assumed to be unique). In this case the variation of h to a new value h' means that either we admit 
that the temperature is actually slightly non-zero, and the system passes to a new ground state configuration 
(corresponding to h') via thermally activated flips, or at each variation of h the ground state configuration is 
achieved by cooling down the system (at fixed h') from the high-temperature disordered state. 

The alternative procedure of the variation of h is essentially different. Here one has to fix the starting 
value of h (e.g. h = 0), and then at any variation of h the spins are just following the directions of their 
local fields (which consist of the contributions due to the regular ferromagnetic interactions with the nearest 
neighbors plus the local value of the random field). In this procedure all thermal "jumps over energy barriers" 
are prohibited, and therefore, most probably, due to the variation of h the system will get stuck in one of the 
metastable states. 

Of course, from the experimental point of view none of these two procedures can be implemented in the 
pure form. The first one (qualitatively, it corresponds to the field cooled (FC) measurements) requires very 
long (formally infinite) waiting time to reach the thermal equilibrium at any variation of h, while the second 
one (it corresponds to the zero field cooled (ZFC) measurements) requires to suppress all thermally activated 
processes which can be achieved only if the temperature is exactly equal to zero. 

In the next Section using simple physical arguments it will be shown that if the strength of the random 
fields is not too small then there exists a region of parameters where the critical fluctuations arc irrelevant, and 
therefore the thermodynamic properties of the system can be described in terms of the saddle-point equations 
of the corresponding (random) Ginzburg-Landau Hamiltonian. It is in this region that the ferromagnetic 
order breaks down. Moreover, using simple scaling arguments it will be shown that in this case the finite 
temperature phase transition in this system is equivalent to the low-temperature phase transition at variation 
of the strength of the random fields. 

In Section 3 the detailed study of the zero-temperature order-disorder phase transition will be done 
in terms of the discrete Ising model on the three-dimensional cubic lattice. Here we propose a modified 
(slightly improved) version of the mean-field approach which takes into account nearest-neighbors two-spins 
correlations. Then, considering the two regimes: increasing and decreasing strength of the random fields, and 
supposing that all thermally activated spin flips are suppressed, it will be demonstrated that the magnetization 
curve m(h) reveals the hysteresis loop. In particular, it is characterized by two different critical points 
hi ' < hi + \ where hi ' is the strength of the random fields at which m(h) becomes non-zero at decreasing 
h, and hi^ is the one at which m{h) becomes zero at increasing h. 

2 Heuristic arguments 

For the sake of generality, before coming down to the dimension three, first, let us consider the continuous 
version of the D-dimcnsional Ising model in terms of the scalar field Ginzburg-Landau Hamiltonian: 



Here r = (T c — T)/T c (r <C 1) is the reduced temperature parameter (for simplicity in what follows it will be 
supposed that T c = 1). According to this definition, positive values of r correspond to the low-temperature 
ferromagnetic state. Random fields h(x) are described by the symmetric Gaussian distribution, 



where h is the parameter which describes the strength of the random field, and po is an irrelevant normalization 
constant. For a given realization of the random fields the partition function of this system is obtained by the 
integration over all configurations of the scalar fields </>(x) 




(2) 




(3) 




(4) 
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It is well known that if we consider such system at temperatures not too close to the critical point (so that the 
parameter r is not too small) the thermal critical fluctuations are irrelevant, and the leading contributions 
to the partition function, eq.(4), are coming from the minima of the Hamiltonian, eq.(2), described by the 
saddle-point equation 



-A^x)-r0(x)+# 3 (x) = /i(x) 



(5) 



In the pure system (at h(x) = 0) the restriction on the value of the parameter r is given by the Ginzburg- 
Landau condition 



r » g 2 /( 4 - D ) = 



tgl 



(6) 



It is clear that this requirement makes sense only if the coupling parameter g is sufficiently small. Of course, 
formally for the three-dimensional Hamiltonian, cq.(2), considered as the continuous limit representation of 
the original Ising model this is not true (unlike the corresponding systems in dimensions close to four, where 
the effective (renormalized) value of the coupling parameter g ~ e = (4 — D) <C 1). Nevertheless, here we 
are going to consider the random field Ginzburg-Landau theory, eq.(2) in which the coupling parameter g 
is assumed to be sufficiently small, hoping that at the qualitative level the behavior of the system is not so 
much sensitive to the actual value of this parameter. 

The only relevant spatial scale in the system described by the Hamiltonian (2) is the correlation length 
R c (t) which under condition (6) has the scaling 



Rc(T) 



-1/2 



(7) 



In the absence of the random fields the ferromagnetic ground state of the system is characterized by the homo- 
geneous configuration </> = \frfg. To study the effects produced by the random fields on the ferromagnetic 
state of the system let us perform the following space and fields rescaling 



0(x) 



1/2 



(f(x/R c ), 



(8) 



In terms of the new fields ip(z), where z = x/i? c = r^x, the system is described by the rescaled Hamiltonian 



(4-P) 



H 



ip(z),h(z) 
where the rescaled random fields 

h(x/R c ) 

arc described by the distribution function 



\ (Vv>(z)) 2 - ±y> 2 (z) + i^ 4 (z) - \/^/>(z) y -(z) 



1 



<i\<Rc 



d°x' /i(x') 



P[h( X )]=p oeX p(-^ m J d D zh 2 (z)^ 



(9) 



(10) 



(11) 



(note that the integration here involves the ultraviolet cutoff lengthscale equal to one). Redefining the random 
fields again 



h(z)= , -£(z) 
instead of eq.(4) we obtain the partition function 

Z[Z(z)] - |^(z)exp(-/3^b(zU(z)] 



(12) 



(13) 



which is controlled by the effective "inverse temperature" 
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(4-Q) 

3,1 <14) 



(tgl is the Ginzburg-Landau temperature, eq(6)) and which is defined by the new effective Hamiltonian 



H[<p(z),£(z)] = J d D z 



\ (V^(z)) 2 - i^ 2 (z) + \p\z) - C(zMz) 



(15) 



which contains no parameters. The new random fields £(z) are described by the Gaussian distribution 



P[£(z)] = Po cxp^-^ | ^£ 2 (z)^) , (16) 



characterized by the mean square value 

A 2 _ (ff^) fl7 N 
A - T (6-D)/2 W 

The ground state configurations in terms of the new fields <p(z) are defined by the saddle-point equation 

-A^(z)-^(z)+^ 3 (z) = £(z) (18) 

In terms of the new fields, <p(z) and £(z), the transition from the ordered (ferromagnetic) to the disordered 
(paramagnetic) state looks as follows. In the absence of the random fields (£ = 0) the ground state is given 
by the trivial ferromagnetic solution (fo(z) = 1 (or <po(z) = —1). The presence of weak random fields (at 
A< 1) introduces only small perturbations to this solution. However, if we increase the effective strength 
of £(z), which is controlled by the parameter A, eq.(17), the ferromagnetic configuration is getting more and 
more perturbed, and finally, at at a certain critical value A c the ferromagnetic ordering is destroyed. 

One can note two important points in the above scenario: 

First. According eq.(17), for fixed value of the parameter h (which is the strength of the original random 
fields /i(x)), the effective strength of the random fields £(z) is controlled by the temperature parameter r. 
According to its definition (r = (1 — T)), increasing the temperature T of the system means decreasing r, 
which in turn, produces increasing of the parameter A In other words, variations in temperature turn into 
variations of the effective strength of the random fields. 

Second. According to the saddle-point equation (18) the critical value of A c is of the order of one. This 
means that the transition takes place at 

t ~ t c = (gh 2 )V^l (19) 
At this point the value of the effective "inverse temperature" parameter (3, eq.(14) is 

(4--D) 

Thus, if the strength of random fields is not too small: 

h » K{g) - gT^l (21) 
at dimensions D < 4 in the vicinity of the phase transition (when r <~ r c ) we have 

2(4-D) 

/ h \"(5 r W 

" ~ * = im) 1 (22> 

In other words, in terms of the representation, eqs.(13)-(17), the original ferromagnetic-paramagnetic phase 
transition at any point of the critical curve h c (T) on the left of the point (/i», T») (Figure 1) is equivalent to 
the low-temperature order-disorder transition (far left extreme of the curve h c (T)) at variation of the strength 
of the random fields. In particular, in the dimension D = 3, 
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tgl = g 2 



* = (t) 



2/3 



(23) 



/i*(flO 



fJ 



The conditions r c ^> tql, Pc ^> 1 and h ^ discussed above are automatically satisfied by the only 
restriction on the strength of the random fields, 



h > g 



(24) 



In the next section to understand the nature of the low temperature phase transition, we are going to consider 
its extreme version in the zero temperature limit. In this case it is natural to consider the original discrete 
Ising model on a lattice instead of its continuous limit representation. 




Figure 1: Qualitative phase diagram of the random field Ising model 



3 Zero-temperature phase transition 



Let us consider the three-dimensional random field ferromagnetic Ising model on the cubic lattice with the 
Hamiltonian 



- N N 



12 

<?'#j> * 

Here < j ^ j > denotes the pairs of nearest neighbors, and the random fields hi are described by the 
independent Gaussian distributions 

N 

W h [h 1 ,h 2 ,...,h N ] = J] V h (hi) (26) 

i=l 

where 

1 cx (_g 



V h (hi) = -j— exp ( ) (27) 
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The ground state spin configuration of this system is defined by the set of N the conditions 



Sign - 



a«=l 



o-q, + h, 



) 



(28) 



where for every site i the summation goes over its 6 nearest neighbors. 

First, let us discuss at the qualitative level what is going on if we increase the parameter h from zero to 
large values. At h = the state of the system is ferromagnetic: the solution of the above equations (28) is 
trivial, <7j = +1 (of course, there is another equivalent solution <7j = — 1, but in what follows we will suppose 
that in the ferromagnetic state the spins are directed "up"). At non-zero h <C 1 almost all the spins are 
directed "up", but in rare cases, when at a given site i the value of the local field is sufficiently negative, 
hi < — 1, the i-th spin is flipped "down". According to eq.(27), the probability of this event (which is equal 
to the concentration of the negative spins) is 



is the complementary error function (which exponentially goes to zero at z ^> 1). 

Note that in the above consideration we have neglected the configurations in which the flipped down 
spins appear to be the nearest neighbors. It is clear that at small h the concentration p(h) is exponentially 
small and therefore the probability of these events is negligible. However, when increasing the value of the 
parameter h the probability to find two or more neighboring "down" spins becomes non-small, and it is due 
to these configuration that the situation becomes rather sophisticated. First of all, one can note that if at 
neighboring sites the values of the random fields appear to be non small (and negative), then the stable 
spin configurations can become ambiguous. As an example let us consider two neighboring sites i and i + 1 
(surrounded by the spins which all are directed "up") with the values of the random fields —1 < hi < —2/3 
and —1 < h i+ i < —2/3. One can easily see that in this case both the configuration Oi = <r i+ i = +1 and the 
configuration Ui — <r i+ i = — 1 satisfy the stability conditions eq.(28). It is obvious that one could find similar 
phenomena in clusters consisting of larger number of spins. 

The effects of the ambiguity of the stationary spin configurations due to the presence of the random fields 
and their (non-perturbative) contributions to the finite-temperature thermodynamics have been discussed 
earlier in terms of the continuous Ginzburg-Landau representation of the RFIM [22] . Here we will describe 
the consequences of the ambiguity of the stationary states for the phase transition at the zero temperatures. 
At the qualitative level it is tempting to suggest that the presence of such kind of phenomena could be the 
origin of the phase transition of the first order. Indeed, if at increasing the strength of the random field the 
system gets stuck at ferromagnetic configurations, an alternative disordered states could become energetically 
preferable before the ferromagnetic state would become unstable. Note that unlike the second-order phase 
transitions which are characterized by the divergence of the correlation length (and which require to find the 
way to study the system at large scales) the first-order phase transition is characterized by a finite value of 
the correlation length. For that reason one can hope that taking into account correlations only of the order 
of the lattice spacing could still give reasonable results. 

3.1 Mean-field approach 

First, as a matter of simple "warming up exercise", let us consider how the the zero-temperature phase 
transition looks like if we study it in terms of the ordinary mean-field approach (which neglects the spin-spin 
correlations at all). 

Let us denote by x the probability that a given spin i points " up" . Then the probability to be " down" is 
or course equal to (1 — x). In this case the ferromagnetic order parameter which is the global magnetization 
of the system is 




(29) 



where 




dt exp (— t 2 ) 



(30) 




(31) 
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Figure 2: Mean field solution for the order ferromagnetic order parameter m(h) (h c ~ 0.696) 



A given spin is surrounded by the six nearest neighbors. If all the neighbors are pointing "up" (the probability 
of this configuration is equal to x 6 ) then Oi — +1 provided the local field hi > — 1. If one of the neighbors 
is "down" (the probability is 6x 5 (l — x)), then <Ji = +1 provided the local field hi > —2/3. If two of the 
neighbors are "down" (the probability is A^x 4 (l — x) 2 ), then <7j = +1 provided the local field hi > —1/3. 
Etc. Collecting all these configurations one easily gets the following self-consistent equation for x: 



x = x 6 ■ P 6 (h) + 6x 5 (l - x) ■ P 5 (h) + 15x 4 (l - x) 2 ■ Pi{h) + 20x 3 (l - xf ■ P 3 (h) + 
+ I5x 2 {l - a;) 4 ■ P 2 (h) + 6x(l - x) 5 ■ P^h) + (1 - a;) 6 ■ P (h) 

Here we have introduced the notations 



(32) 



Pe(h) 
P 5 (h) 
P4(h) 
Ps(h) 
Pz(h) 
P 1 (h) 
Po(h) 



dyV h (y) = 



2/3 



dyV h (y) 



1/3 

dyPh(y) 



1/3 



dyV h {y) 



2/3 



dyV h {y) 



2 \y/2h 



1 



dyVh(y) = 



-Erfc 



1 - -Erfc 
2 



3V2h 
1 

3\/2/i 



dyV h (y) = 



^Erfc 



-Erfc 



^Erfc 



1 



3V2h 
2 

3V2h 
1 

V2h 



(33) 



The solution of eq.(32) for the ferromagnetic order parameter is shown in Figure 2. We see that the phase 
transition from the ferromagnetic (m ^ 0) to the paramagnetic (m = 0) state takes place at h c ~ 0.696. One 
can easily check that in the vicinity of the critical point, as h — > h c , the the ferromagnetic order parameter 
vanishes according to the the scaling law: 



m(h) - {hc-hf 



(34) 



described by the usual mean-field critical exponent [3 — 0.5 (see Figure 3). This, of course, is not surprising 
because if in eq.(32) we substitute x = (m+ l)/2, the resulting equation for m (whatever complicated it may 
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Figure 3: Log-log plot for the ferromagnetic order parameter in the vicinity of the critical point 



look like) will be still the mean field equation for the order parameter, and its development in powers of small 
m(h) in the vicinity of the critical point can not give anything else but mean- field critical exponent. 

It is clear that the consideration presented above can not pretend to describe the true phase transition 
which takes place in the system under consideration. The reason for that is obvious: in such king of the 
mean-field calculations one misses all spin-spin correlations, which, as we have discussed in the beginning of 
this section, are crucial for the statistical properties of the present system. 

Below we are going to modify the above approach so that it would take into account correlations between 
pairs of nearest neighbors spins. It turns out that even such rather limited improvement is sufficient to 
intoduce rather dramatic changes in the scenario of the phase transition. 



o 



Figure 4: Two neighboring sites are surrounded by (5+5) nearest neighbors 



3.2 Modified mean-field approach 

Let us consider a pair of nearest neighbor spins o\ and 02 and let us denote by p(Tt), p(T4), p(4T) an( i p(ll) 
the probabilities of the (01, 02) configurations (+1,+1), (+1,-1), (— 1,+1) and (—1,-1) correspondingly. 
These four probabilities are, of course, bounded by the condition p(TT) +P(t-L) +p(lt) +p(-U) = 1- As before, 
by x we denote the probability for a given spin to be "up". One can easily see that 

x = p(TT) + p(U) (35) 
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The two spins o~i and 02 are surrounded by (5+5) neighbors (Figure 4), and their orientations are defined by 
the equations: 



0"! 



02 = Sign 



Sign \ \ J~] a ai + cr 2 I + hi 
\ f Yl <r aa + <ri) + h 2 



(36) 



where <J ai 2 denote the neighbors of the spins oi^ 2 . 

The idea is to compute p(TT) an d P(T-l) by summing over all possible configurations of these neighboring 
spins (with the corresponding probabilities defined by the parameter x) and integrating over local fields hi 
and h 2 (with the probability distribution Vh(hi)Vh(h 2 )). Then, substituting p(TT) an d p(tl) (which will be 
the functions of x and h) into eq.(35), we will get a self-consistent equation for x. It turns out, however, that 
this program can not be implemented directly, because for every configuration of surrounding spins there are 
finite regions in the plane [hi, hi), where the orientations of spins o~i and 02 are not uniquely defined by the 
given values of hi and hi. 




o c 
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Figure 5: The structure diagram of the orientations of two neighboring spins for the case when all their (5+5) 
nearest neighbors are directed "up" 



As an example let us consider again the simplest situation when all ten neighboring spins are directed " up" 
(which has the probability a; 10 ). In this case the condition to have both spins" up" is hi > — 1 and hi > — 1 (so 
that the corresponding contribution to p(TT) must be given by the integration oi'Ph(hi)'Ph(h 2 ) in the sector 
(hi > —1, h 2 > —1))- The contribution to the probability p(H) is given by the sector (hi > —2/3, h 2 < —1). 
The contribution to the probability p(ll) is given by the sector (hi < —1, hi > —2/3). And finally the 
contribution to the probability p([ [) is given by the sector (hi < —2/3, h 2 < —2/3), which overlaps with the 
sector corresponding to p(TT)- Thus in the square (—1 < hi < —2/3; — 1 < h 2 < —2/3) we are facing a kind 
of frustration (Figure 5). 

One can propose two ways out of this situation. First of all, one can, of course, compare the energies of the 
two configurations, and then, if the spin orientations always correspond to the ground state, the ambiguity 
will be lifted (the configuration with both spins "up" has lower energy in the triangle above the dotted line 
inside the "frustrated square" in Figure 5, and the one with both spins "down" has lower energy below the 
dotted line). However, if we keep in mind the scenario in which in the course of variations of the external 
parameter h each spin just follow the direction of the local field, while all thermally activated "jumps over 
barriers" are suppressed, then we have to conclude that inside the "frustrated" region the orientation of 
the two spins must depend on the history. Namely, for every point in the plane (hi,h 2 ) we can assign the 
"trajectory" which would demonstrate where this particular values of the random fields came from. For 
instance, in the situation when we study the evolution of the system for increasing value of h, (starting from 
zero) a particular point (hi, hi) has the "trajectory" which is the straight line connecting this point with the 
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Figure 6: The orientation diagram of two neighboring spins for the case when all their (5+5) nearest neighbors 
are directed " up" : (a) at increasing and (b) at decreasing strength h of the random fields 



origin. In the opposite case of decreasing h (starting from infinity), the "trajectory" of the point (Aij^a) is 
the straight line which comes from infinity and directed towards the origin. 

Following these rules of the game the ambiguity of the spin orientations is lifted provided one fixes which 
process is actually under study: increasing or decreasing of h. In the study of the ferromagnetic — ► paramag- 
netic transition (increasing h) one can easily note that the "trajectory" to any point inside the "frustrated" 
square (—1 < hi < —2/3; — 1 < h 2 < —2/3) comes from the region where both spins, o\ and 02, are directed 
"up". Thus, in this case the square "belong" to the probability p(TT) (Figure 6(a)). In the reversed process 
of decreasing h (the paramagnetic — > ferromagnetic transition) the only way to arrive inside this square is 
to come from the region where both spins are " down" , and in this case the square belong to the probability 
p{ll) (Figure 6(b)). 

For better understanding how it works let us consider other examples. In the case when all 10 surrounding 
spins are directed " down" , the " orientation diagram" for the spins o\ and 01 is represented in Figure 7. Here 
the "frustrated square" is located in the region (2/3 < h\ < 1; 2/3 < hi < 1). Now, at increasing h one 
arrives into this square from the region where both spins are " down" and therefore in this case it belongs to 
the probability p(H) (Figure 7(a)). Apparently, in the reversed process of decreasing h the square belong to 
the probability p( ff ) (Figure 7(b)). 

Let us consider the configuration in which all 5 neighbors of the spin a± are " up" , while all 5 neighbors of 
the spin 02 are "down". Using eqs.(36) one can easily build the corresponding orientation diagram (Figure 
8). We see that in this case both at increasing and decreasing h one can arrive into the left-down triangle of 
the "frustrated square" only by passing the region where both spins are "down" . Similarly, both at increasing 
and decreasing h one can arrive into the right-up triangle only by passing the region where both spins are 
"up". 

Finally, let us consider the configuration in which all 5 neighbors of the spin a\ are " up" , and only 4 
neighbors of the spin 02 are "down". Here the division of the "frustrated square" between p(TT) and p(H) 
is somewhat more complicated (Figure 9). At increasing h one arrives into the section of the square below 
the line h2{h\) = —\h\ by passing through the region where both spins are "down" and one arrives into the 
section above this line by passing through the region where both spins are "up" (Figure 9(a)). On the other 
hand, at decreasing h one arrives into the section of the square below the line /12C11) = — § h% by passing 
through the region where both spins are "down" and one arrives into the section above this line by passing 
through the region where both spins are "up" (Figure 9(b)). 

Now, with some patience and perseverance one can construct the structure of the frustrated squares for 
all 36 "analytically different" (in terms of eqs.(36)) configurations of the surrounding spins. The resulting 
orientation diagrams for increasing and decreasing h are shown in Figures 10 and 11 correspondingly. 
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Figure 7: The orientation diagram of two neighboring spins for the case when all their (5+5) nearest neighbors 
are directed "down": (a) at increasing and (b) at decreasing strength h of the random fields 



Finally, after the regions corresponding to "up-up" and "up-down" orientations of the spins G\ and a 2 
are unambiguously defined, the derivation of the equations for the probabilities p( J \'\)(x 1 h) and p{][){x, h) is 
straitforward (although slightly cumbersome). The explicit forms of the equation (35) for increasing and for 
decreasing h is given in the Appendix. 

The solution of these equations in terms of the ferromagnetic order parameter m(h) is represented in 
Figure 12. We see that the behavior of the magnetization exhibits clear hysteresis phenomenon. 



4 Discussion 

In this paper the zero-temperature phase transition has been studied under the assumption that at any 
variations of the strength h of the random fields the transformations of the spin configurations go under 
constrain that all thermally activated spin flips are suppressed. It should be stressed that this situation is 
essentially different from the true equilibrium phase transition where at any given h the system is supposed to 
be in the ground state (and which, in my view, is experimentally inaccessible, et least in the low-temperature 
limit). 

Comparing two possible scenarios of the order-disorder phase transition: continuous (the second order) 
and discontinuous (the first order), one should note that unlike the continuous transition characterized by 
the divergence of the correlation length, at the first-order phase transition the correlation length remains 
finite. Thus, admitting that the transition is discontinuous one can hope that a theory which takes into 
account spin-spin correlations only at a limited scale (e.g. of the order of the lattice spacing) still would give 
qualitatively correct description of the phase transition. 

Usual mean-field theory is useless here because it does not take into account spin-spin correlations at all. 
On the other hand, it has been demonstarted in this paper that even rather limited improvement of the mean- 
field approach, which takes into account two nearest neighbor spins correlations, produces rather dramatic 
effect on the scenario of the phase transition. After this modification the value of the ferromagnetic order 
parameter m(h) as the function of the strength of the random fields becomes history dependent exhibiting 
clear hysteresis phenomenon (Figure 12). 
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Figure 8: The orientation diagram of two neighboring spins for the case when all 5 neighbors of the spin <j\ 
are directed " up" , while all 5 neighbors of the spin 02 are directed " down" 




(a) <b) 



Figure 9: The orientation diagram of two neighboring spins for the case when all 5 neighbors of the spin u\ 
are directed " up" , and 4 neighbors of the spin 02 are directed " down" : (a) at increasing and (b) at decreasing 
strength h of the random fields. The lines are: [1] h2{h\) — — (l/2)/ii; [2] h2(h\) — — (2/3)/ii 
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Figure 10: The structure of the "frustrated squares" at increasing strength h of the random fields. The 
notation (m, n) inside the square indicate that it corresponds to the configuration with m neighbors of the 
spin cti and n neighbors of the spin er 2 are directed "up". Rose regions corresponds to the state with both 
spins <7i and cr 2 are directed "up". Green regions are the ones where both spins o\ and 172 are directed 
"down". The lines are: [1] h 2 (hi) = -(1/2)%; [2] h 2 (hi) = -hi; [3] h 2 (hi) = -2hi. 




Figure 11: The structure of the "frustrated squares" at decreasing strength h of the random fields. Notations 
are the same as in Figure 10. The lines are: [1] h2(hi) — —(l/3)hi; [2] ^2(^1) = — (l/2)fti; [3] /i 2 (^i) = 
-(2/3)fti; [4] hi{h!) = -hi; [5] h 2 (hi) - -(3/2)%; [6] h 2 (hi) = -2%; [7] h 2 (hi) = -3% 
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Figure 12: Ferromagnetic order parameter m(h) at increasing (red line) and decreasing (green line) strength 
h of the random fields. ~ 0.722; h ( ~ ] ~ 0.622. 



Unfortunately in the framework of the present theory it is still difficult to to make a definite conclusion 
about the nature of the phase transition. The existence of the hysteresis loop (together with the clear 
understanding of the physical mechanism in its origin) is the strong argument in favor of the first-order phase 
transition. On the other hand, each of the curves m + (h) (at increasing h) and m_(/i) (at decreasing h) 
in the Figure 12 demonstrates the continuous transition which, of course, makes no sense, because at such 
transition the correlation length diverges, while the present theory takes into account the correlations only 
of the order of the lattice spacing. In this sense the present theory is not self-consistent, and it would be 
reasonable to expect that the presence of the continuous transitions at he is not more than just an artifact 
of the proposed approach. These issues require further detailed studies. 



APPENDIX 

Let us denote by p+(tt) an( l P-(TT) t ne probabilities for the two spins a\ and a 2 to be both "up" at 
increasing and at decreasing variations of h correspondingly. Note that the probability p(TI) is the same in 
both (increasing and decreasing) cases. Then, using the orientation diagrams shown in Figures 10 and 11 
(after some work) we get: 
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p+(tf)(x,h) = P 6 (ft)P 6 (ft)z 10 + 10P 6 (ft)P 5 (ft)z 9 (l - xf 

+ [20P 6 (ft)P 4 (ft) + 25P 5 (ft)P 5 (ft)] x s {l - xf + [20P 6 (ft)P 3 (ft) + 100P 5 (ft)P 4 (ft)] x 7 (l - xf 
P 6 {h)P 2 (h) - D { +\h)\ + 100P 5 (ft)P 3 (ft) + 100P 4 (ft)P 4 (ft) 



+ 



+ 



10 



x 6 (l-xf 



2 [P 6 (ft)Pi(ft) - 0.5(P 6 (ft) - P 5 (ft))(Pi(ft) - Po(h))} + 50 [P 5 (ft)P 2 (ft) - 0.5(P 5 (ft) - P 4 (ft))(P 2 (ft) - Pi (ft))] 



+200 [P 4 (ft)P 3 (ft) - 0.5(P 4 (ft) - P 3 (h))(P 3 (h) - P 2 (h))} 



x 5 (l-xf 



+ 



10 



P 5 (h)Pi(h) - (P 5 (h) - P 4 (/i))(Pi(/i) - Po(h)) + D ( +\h) + 100 [P 4 (ft)P 2 (ft) - (P 4 (h) - P 3 (ft))(P 2 (ft) - P(ft))] 



+ 100 [P 3 (h)P 3 (h) - (P 3 (h) - P 2 (h))(P 3 (h) - P 2 (h))} 



x 4 {l-xf 



20 [P 4 (ft)Pi(ft) - (P 4 (ft) - P 3 (ft))(Pi(ft) - P (ft))] + 100 [P 3 (ft)P 2 (ft) - (P 3 (ft) - P 2 (h))(P 2 (h) - P^h))} 



+ 20 (P 3 (/i)Pi(/i) - (P 3 (ft) - P 2 (ft))(Pi(ft) - P (ft))) + 25(P 2 (ft)P 2 (ft) - (P2W - Pi(h))(P 2 (h) - P(ft))) 



x 3 (l-x) 7 



x 2 (l-2;) 8 



+ 10 [P 2 (ft)Pi(ft) - (P 2 (ft) - Pi(fc))(Pi(/i) - P (ft))] ^(1 - a;) 9 
+ [Pi(ft)P(ft) - (Pi (ft) - P (ft))(Pi(ft) - P (ft))] (1 - x) 10 



(37) 



p-(TT)(^,/i) = 

= P 6 (h)p e (h) - (P 6 (ft) - P 5 (ft))(P 6 (ft) - P 5 (h)) 



x w + 10 



P 6 (ft)P 5 (ft) - (P 6 (h) - P 5 (ft))(P 5 (ft) - P 4 (ft)) 



+ 



20 [P 6 (ft)P 4 (ft) - (P e (h) - P 5 (ft))(P 4 (ft) - P 3 (h))} + 25 [P 5 (h)P 5 (h) - (P 5 (ft) - P 4 (h))(P 5 (h) - P 4 (ft))] 



x 9 (i-xy 

x s (l-xf 



+ 



20 



10 



P 6 (h)P 3 (h) - (P 6 (ft) - P 5 (ft))(P 3 (ft) - P 2 {h)) + D ( -\h) 



+ 100 [P 5 (ft)P 4 (ft) - (P 5 (ft) - P 4 (ft))(P 4 (ft) - P 3 (ft))] x 7 (l - xf 



P 6 (h)P 2 (h) - (P 6 (ft) - P 5 (ft))(P 2 (ft) - P 1 (h)) + Dfc\h) 



+ 100 



P 5 (ft)P 3 (ft) - (P 5 (h) - P 4 (h))(P 3 (h) - P 2 {h)) + D^\h) 



+ 100 [P 4 (h)P 4 (h) - (P 4 (ft) - P 3 (h))(P 4 (h) - P 3 (h))} 



x 6 (l-xf 



+ 



2 [P 6 (ft)Pi(ft) - 0.5(P 6 (ft) - P 5 (ft))(Pi(ft) - P (ft))] + 50 [P 5 (ft)P 2 (ft) - 0.5(P 5 (ft) - P 4 (ft))(P 2 (ft) - Pi (ft))] 



+200 [P 4 (ft)P 3 (ft) - 0.5(P 4 (ft) - P 3 (ft))(P 3 (ft) - P 2 (ft))] 



x 5 (l -x) 5 



10 



20 



P 5 (h)P 1 (h)-D^\h) 



100 



P 4 (ft)P 2 (ft) - /^(ft) + 100P 3 (ft)P 3 (ft) 



a; 4 (l-a;) 6 



P 4 (ft)Pi(ft) - ^ 2 '(ft) + 100P 3 (ft)P 2 (ft) 



x 3 (l - a;) 7 



20P 3 (ft)Pi(ft) + 25P 2 (ft)P 2 (ft)] x 2 (l - x) 8 + 10P 2 (ft)Pi(ft)a; 1 (l - ^) 9 + Pi(ft)Pi(ft)(l - x) 



10 



(38) 
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p(U)(x,h) = P 5 (h)P (h)x w + [P 5 (/i)Pi(/i) + 5P 4 {h)P (h)]x 9 (l - xf 

+ [WP 5 (h)P 2 (h) + 25P 4 (h)P 1 (h) + 10P 3 (h)P Q (h)} .t 8 (1 - xf 

+ [WP 5 {h)P 3 (h) + 50P 4 (h)P 2 (h) + 50P 3 (/i)Pi(ft) + 10P 2 (h)P (h)} x 7 (l - xf 

+ [5P 5 {h)P 4 (h) + 50P 4 (h)P 3 {h) + 100P 3 (h)P 2 {h) + 50P 2 (/i)Pi(/i) + 5Pi(h)P (h)] x 6 {l - xf 

+ [P 5 (h)P b {h) + 25P i (h)P i {h) + l00P 3 (h)P 3 {h) + li){)P 2 (h)P 2 {h) + 25P 1 {h)P 1 (h) + P (h)P {h)] x 5 (l - xf 

+ [5P 4 (h)P 5 (h) + 50P 3 (h)P 4 (h) + W0P 2 (h)P 3 {h) + 50Pi(/i)P 2 (/i) + 5P (/i)Pi(/i)] x 4 (l - xf 

+ [WP 3 {h)P 5 (h) + 50P 2 {h)P 4 (h) + 50Pi(/i)P 3 (/i) + 10P {h)P 2 (h)] x 3 (l - xf 

+ [li)P 2 {h)P 5 (h) + 25P 1 {h)P 4 (h) + 10P (h)P 3 {h)} x 2 (l - xf 

+ [5P!(/i)P 5 (/i) + 5P (h)P 4 (h)} x\l - xf + P (h)P 5 (h)(l - xf° 

(39) 

where the functions Pk(h) (k = 0, 6) are defined in cq.(33) and 

r-2/3 ,- yi /2 

- / rfyi^(yi) / d^G/a) 

J-l Jl/3 
pi /--1/3 

J2/3 J-V2/2 
r-2/3 ,1/3 

^raW = / d yi P h ( yi ) / dy 2 V h (y 2 ) 
J-l J-yi/3 

2^3 2/3 

= / dyin(wi) / ^2^(2/2 ) 

J-l J-2 yi /3 
r-1/3 ,1/3 

-D42 C 1 ) = / dyiPhivi) / dy 2 V h (y 2 ) (40) 
J-2/3 J-yi/2 

The ferromagnetic order parameters m±(h) as the function of ft, are obtained from the relation 

m±(/i) = 2x±(ft) - 1 (41) 
where and X-(h) are the corresponding solutions of the equations 

x = P±(tf)(x,h) + p(U)(x,h) (42) 
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